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To discuss the salient role of the statistical memory effects in the human brain functioning we 
have analyzed a set of stochastic memory quantifiers that reflects the dynamical characteristics of 
neuromagnetic brain responses to a flickering stimulus of different color combinations from a group 
of control subjects which is contrasted with those from a patient with photosensitive epilepsy (PSE). 
We have discovered the emergence of strong memory and the accompanying transition to a regular 
and robust regime of chaotic behavior of the signals in the separate areas for a patient with PSE. 
This finding most likely identifies the regions of the location the protective mechanism in a human 
organism against occurrence of PSE. 

PACS numbers: 05. 45. Tp; 87. 19. La; 89. 75. -k 

Increasing attention has been paid recently to the study of statistical memory effects in random processes that 
originate from nature by means of nonequilibrium statistical physics. The role of memory has its roots in natural 
sciences since 1906 when the famous Russian mathematician Markov wrote his first paper on the theory of Markov 
Random Processes (MRP) His theory is based on the notion of an instant loss of memory from the prehistory 
(mcmoryless property) of random processes. In contrast, there are an abundance of physical phenomena and processes 



which can be characterized by statistical memory effects: kinetic and relaxation processes in gases [2| and plasma 

condensed matter physics (liquids solids a], and superconductivity Q|), astrophysics nuclear physics 
quantum [9] and classical Q physics, to name only a few. At present, we can make use of a variety of statistical 
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methods for the analysis of the memory effects in diverse physical systems. Typical such schemes are Zwanzig-Mori's 
kinetic equations lll| , generalized master equations and corresponding statistical quantifiers 
relation method [13| , the generalized Langevin equation (GLE) 14 1 , etc. 



12j . Lee's recurrence 



In this paper we shall demonstrate that the presence of statistical memory effects is of salient importance for the 
functioning of healthy physiological systems. Particularly, it can imply that the presence of large memory times scales 
in the stochastic dynamics of discrete time series can characterize pathological (or catastrophical) violation of salutary 
dynamic states of the human brain. As an example, we will demonstrate here that the emergence of strong memory 
time scales in the chaotic behavior of neuromagnetic responses of human brain as recorded by MEG is accompanied 
by the likely initiation and the existence of PSE. 

First consider a simplified version of the Markov processes. Let us introduce the conditional probability 
Ki(xi, ti\x2, t2) that x is found in the range {x2,X2 + dx%) at X2, if x had the value x\ at t\. For the Markov 
random process the conditional probability that x lies in the range (x n ,x n + dx n ) at t n given that x had the values 
xi, X2, ■■■x n -i at times t\, £2, ■■■■t n -i depends only on x n ^\ is as follows: K n -i(x\,ti;X2,t2; ...x n -i,t n -i\x n ,t n ) = 
Ki(x n -i,t n -i\x n ,t n ). The equation states that, given the state of a Markov process at some times t Tl -i < t n , the 
forthcoming (future) state of the process at t n is independent of all previous states at prior times. The equation is 
the standard definition of the Markov random process. So, from the physical point of view the Markov process is the 
process without aftereffect. It means that the "future " and the "past" of a process not depend each from other at 
known " present" . 

Measures for memory. One of the first measure of 'memory' in physiological time series that has been studied in 
electroencephalographic (EEG) and magnetoencephalographic (MEG) signals, both of healthy subjects and patients 
(including epilepsy patients) [3] was the detrended- fluctuation analysis (DFA) 

For the quantitative description of statistical memory effects of random processes in the physiological data the use 
of Zwanzi-Mori kinetic equations provides an appropriate and most convenient methodology. In particular, using 
the reasoning put forward in Refs. [171 ] one can obtain the chain of interconnected kinetic equations for the discrete 
time correlation function (TCF) a(t) = M (t) — {Sx(t)Sx(0))/{Sx 2 (0)) of the fluctuation Sx(t) = x(t) — (x(t)), where 
x(t) = [x\\ X2', xn) is a random discrete-time process, i.e., Xj = xj(tj), tj=j r, where r is a discretization time-step, 
j = 1, 2, ...N. This zeroth-order function is then related iteratively to higher order memory functions Mi(t), i = 1,2, .... 
In this approach the set discrete memory functions Mi(t),i — 1, 2, ... (MF's) of i th order together with corresponding 
relaxation parameters quantify the memory effects. The full set of MF's includes all peculiarities of the memory effects 
for real complex systems. For the discrete time series the whole set of functions Mi(t) and relaxation parameters can 
be calculated directly from the experimental data [17|. 

17[ provides the adequate tools to study the role of memory 



Following the reasoning put forward with Refs. 



effects in discrete time complex systems dynamics. The characterization of memory is based on a set of dimensionless 
statistical quantifiers which are capable of measuring the strength of memory that is inherent in the complex dynamics. 
A first such measure is is ei(ui) = {^{uj) / /ii+i^)} 1 / 2 whereas the as second set of measure follows as Si(u>) = \M-(ui)/ 
M' i+1 (io)\. Here, ^i(uo) = \Mi{u>)\ 2 denotes a power spectrum of the corresponding memory function Mi(t), Ml(ui) = 



dMfiuj) / dui and Mj(w) is a Fourier transform of the memory function Mj(t). The measures Ei(u>) are suitable for 
the quantification of the memory on a relative scale whereas the second set Si(u>) proves useful for quantifying the 
amplification of relative memory effects occurring on different complexity levels. Both measures provide statistical 
criteria for the comparison of the relaxation time scales and memory time scales of the process under study. For values 
obeying {e, 6} » 1 one can observe a complex dynamics characterized by short-ranged temporal memory scales. In 
the limit these processes assume a 5-like memory with e, S — > oc. When {e, 6} > 1 one deals with a situation with 
moderate memory strength, and the case with both e, 6 ~ 1 typically constitutes a more regular and robust process 
possessing strong memory features. 
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Figure 1: Time dependence of TCF Mo(t)(i = 0) and first two subordinate MF's Mi(t),i = 1,2 for a healthy subject (No. 
6) (blue) and for a patient (red) with PSE for the SQUID's number n = 10, r = 0, 2 ms. The drastic distinctions of Mi(t) 
in a healthy person as compared to a patient with PSE is clearly detectable. They consist in the appearance of significant 
long-ranged oscillations in a healthy subject and the suppression of high frequency noise in a patient with PSE. 



Experimental data for PSE. Next, we can proceed directly to the analysis of the experimental data: MEG signals 
recorded from a group of nine healthy human subjects and in a patient with (PSE) [18 1. PSE is a common type of 
stimulus-induced epilepsy, defined as recurrent convulsions precipitated by visual stimuli, particularly a flickering light. 
The diagnosis of PSE involves finding paroxysmal spikes on an EEG in response to the intermittent light stimulation. 
To elucidate the color-dependency of PS in normal subjects, brain activities subjected to uniform chromatic flickers 
with whole-scalp MEG has been measured in Ref. (further details of the MEG experiment one can find in [lsI ]. 

Nine-right-handled healthy adults (two females, seven males; age range 22-27years) voluntarily participated. Sub- 
jects were screened for photosensitivity and personal or family history of epilepsy. The experimental procedures 
followed the Declaration of Helsinki and were approved by the National Children's Hospital in Japan. All subjects 
gave their informed consent after the aim and potential risk of the experiment were explained. During the recording, 
the subjects sat in the magnetically shielded room and were instructed to observe visual stimuli passively without 
moving their eyes. 

Stimuli were generated by the two video projectors and delivered to the viewing window in the shield room through 
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an optical fiber bundle. Each projector continuously produced a single color stimulus. Liquid crystal shutters were 
located between the optical device and the projectors. By alternative opening one of the shutters for 50 ms, 10 Hz 
(square-wave) chromatic flicker was produced on the viewing distance of 30 cm. Three color combination were used 
: red-green (R/G), blue-green (B/G), and red-blue (R/B). CIE coordinates were x=0.496, y=0.396 for red; x=0.308, 
y=0.522 for green; and x=0.153, y= 0.122 for blue. All color stimuli had a luminance of 1.6 cd/m 2 in otherwise total 
darkness. In a single trial, the stimulus was presented for 2s and followed by an inter-trial interval of 3s, during which 
no visual stimulus was displayed. In a single session, color combination was fixed. 

Neuromagnetic responses were measured with a 122-channel whole-scalp neoromagnetometer (Neuromag-122; Neu- 
romag Ltd. Finland). The neoromag-122 has 61 sensor locations, each containing two originally oriented planner 
gradiometers coupled to dc-SQUID (superconducting quantum interference device) sensors. The two sensors of each 
location measure two orthogonal tangential derivatives of the brain magnetic field component perpendicular to the 
surface of the sensor array. The planner gradiometers measure the strongest magnetic signals directly above local 
cortical currents. From 200 ms prior responses were analog-filtered (bandpass frequency 0.03-100 Hz) and digitized 
at 0.5 kHz. Eye movements and blinks were monitored by measuring an electro-oculogram. Trials with MEG am- 
plitudes > 3000 fT/cm and/or electro-oculogram amplitudes > 150 /zV were automatically rejected from averaging. 
Trials were repeated until > 80 responses were averaged for each color-combination. The averaged MEG signals were 
digitally lowpass-filtered at 40 Hz, and then the DC offset during the baseline (—100 to ms) was removed. At each 
sensor location, the magnetic waveform amplitude was calculated as the vector sum of the orthogonal components. 
Peak amplitude were normalized within each subject with respect to the subject's maximum amplitude. The latency 
range from —100 to —1100 ms was divided with 100 ms bins. Then, the peak amplitudes were calculated by averaging 
all peak amplitudes within each bin. 

Memory analysis for presence of PSE. 

With our set of Figs. 1-5 we present the results of numerical calculations and the analysis of the experimental 
data within the framework of the nonequilibrium statistical approach for stochastic processes in the discrete complex 
systems 17|. In Figs. 1 - 3 we depict the typical data for one concrete healthy subject (No. 6) in comparison with a 
PSE patient for the case of a Red-Blue (RB) combination of the color stimulus. To make the conclusion about the 
role of the statistical memory effects we also show the averaged data for the whole group of nine healthy subjects 
versus the patient with PSE in Figs. 4, 5. 

Fig. 1 depicts the time dependence of the time correlation function (TCF) Mo(t) and the first two MF's Mi(t), 
i = 1,2 for a healthy subject (No. 6) (blue) versus a patient with PSE (red). The TCF Mq(£) displays long-ranged 
oscillations in the healthy and a sharp decay for the patient with PSE. As one can observe from Figs. 2, where the 
power spectra of TCF and MF's are represented, the fractal dependence at order 0; i.e., fio(u>) ~ uj~ a with a = 1, 74 in 
the TCF of the healthy person (blue) now transforms into a group of peaks corresponding to a, (3, 7, S and 9 rhythms 
in frequency behavior of the subordinate quantifiers fii(u>), with i = 1,2,3. The typical picture in the patient with 
PSE (red) consists in (i) the characteristic absence of the fractal dependence for uo(ui), in (ii) the disappearance of 
the well-defined manifestation of physiological e/m rhythms and (iii) in the appearance of a single spike peak at the 
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Figure 2: Power spectra = 0, 1,2 for the MF's in a healthy person (blue) and in the patient with PSE (red) for sensor 

number n=10 in double-log scale. The spectra in the healthy person(No. 6) demonstrate the presence of electromagnetic waves 
on characteristic frequency scales of a, f3, 7, 8 and 6 rhythms (in ^(oj))- One can observe the noticeable peaks of electromagnetic 
excitations in a patient with PSE near 50 Hz and 100 Hz. The similar peaks are present in many other sensors of the human 
cerebral cortex with PSE. The fractal dependence /io(<^) ~ that typifies a healthy person is absent in a patient with PSE. 
This transition plays a crucial role for the emergence of strong memory in a patient with PSE. 
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Figure 3: The frequency dependence of the first three subordinate statistical quantifiers measuring the strength of memory 
£i(u)),i = 1, 2,3 in the healthy person (No. 6) (blue) and for a patient with PSE (red) for SQUID's number n — 10. A distinct 
reduction by a factor 1/56, 5 in the zero-frequency value ei(u> = 0) occurs from the healthy person as compared to the patient 
with PSE. This feature quantifies the emergence of strong memory in a subject with PSE. It is further accompanied by a 
noticeable disappearance of sharp e/m excitations at low frequencies and by the appearance of high frequency noise. 

frequency of 101,5 Hz in the all spectra and for all sensors n. 

The most instructive singularities in the frequency dependence of the first three points of the measure of memory 
Ei(u>), i = l,2,3 (Figs. 3) are as follows. In the healthy person we observe: the fractal dependence in the low frequency 
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Figure 4: The topographic dependence of the information measure for memory ei(u> = 0; n) in the healthy person(Mue) (at fixed 
n, the mean value for the whole group of the 9 control subjects)is compared with the patient with PSE (red), n = 1,2, 3, ...61 is 
the sensor number on the human cerebral cortex. The crucial role of the strong memory for n — 10, 46, 51, 53 and 59 is clearly 
detectable. All sensors depicting ei(0u> = 0; n) surely demonstrate the emergence of statistical memory effects in the chaotic 
behavior of magnetic signals. Nevertheless, the role strong memory effects, i.e. minimum values for E\(uj = 0; n), appreciable 
increases in the patient in the sensors with numbers n = 10, 46, 51, 53 and 59. 

area (uj < 50Hz) £\(lo) ~ cj~' 3 with (3 = 1,67 , the specific behavior with E2(oj = 0) — > and 2 single peaks 

in the area of the frequencies of the brain rhythms for the third point £3(10). This behavior is characteristic only of 
the healthy subjects. The role of increasing memory and the persistent transition from a more random (healthy) into 
a robust, more regular regime of the underlying chaotic process at all three subordinate measures £i(ui),i = 1,2,3 
is clearly detectable in the patient with PSE. The crucial role of the strong memory at the first level, i.e. for e\ is 
reflected by a decrease of the memory measure £\{u)q = 0) by a factor of ca. 56. Moreover, there occurs a drastic 
change of the frequency spectra for Ei(lo) and £3(0;). 

The topographic dependence of £i(oj = 0; n) depicted in Fig. 4 demonstrates the existence of the long-range time 
correlation accompanied by a pronounced increase of the role of the statistical memory effects in all MEG's sensors 
with sensor numbers n — 1,2, ...61 in a patient with PSE as compared to healthy persons. There occurs about one 
magnitude of difference between healthy subject and subject with PSE. 

To specify the role of the strong memory we further study the spatial dependence in terms of a novel information 
measure, the index of memory, which is defined by: 

i/(n) - S^ ealthy (0;n)/SP aUent (Q;n), 

see in Fig. 5. This measure quantifies the sharp revising of memory effects in individual MEG sensors in the patient 
with PSE versus the healthy group. The sharp increase of the role of the memory effects in the stochastic behavior 
of the magnetic signals is clearly visible for sensor numbers n — 10, 46, 51, 53 and 59. The observed points of MEG 
sensors locate the regions of a protective mechanism against PSE in a human organism: frontal (sensor 10), occipital 
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Figure 5: Topographic dependence of the index v(n) (see text) versus sensor n, being the SQUID's number on the human 
cerebral cortex. This indicator measures the amplification of the role of the memory effects. The sharp increase of v(n) for 
n = 10, 46, 51, 53 and 59 characterizes the noticeable increase of the memory effects in the chaotic behavior of magnetic signals 
in the patient with PSE and thus emphasizes the crucial role of the location and the pathological mechanism of PSE. 



(sensors 46, 51 and 53) and right parietal (sensor 59) regions. The early activity in these sensors may reflect the 
protective mechanism that suppresses cortical hyperactivity due to chromatic flickering. 

One might remark that some earlier steps towards the understanding the normal and diseased human brain have 
already been set in other fields of science such as neurology, clinical neurophysiology, neuroscience and so on. The 
numerous studies applying linear and nonlinear time series analysis to EEG and MEG in epileptic patients are 
discussed in details in Refs. [lij ] with taking into account the neurophysiological basis of epilepsy, in particular 
photosensitive epilepsy. Specifically, the results of [18j suggested that a significant nonlinear structure was evident in 
the MEG signals for control subjects, whereas nonlinearity was not detected for the patient. In addition, the couplings 
between distant cortical regions were found to be greater for control subjects. The important role of combinational 
chromatic sensitivity in sustained cortical excitation was also confirmed. These prior finding lead to the hypothesis 
that the healthy human brain is most likely equipped with significantly nonlinear neuronal processing reflecting an 
inherent mechanism defending against hyper-excitation to chromatic flickering stimulus, and such nonlinear mechanism 
is likely to be impaired for a patient with PSE. 

Conclusions. This study of the chaotic behavior of the neuromagnetic signals of a human MEG's with PSE and 
in a group of healthy subjects elucidates the role of the statistical memory as an important criterion, measuring 
the functioning of the human brain. Even an insignificant amplification of the memory effects tests the pathological 
changes in the brain of a patient with PSE. The pronounced sharp increases of memory effects in our set of statistical 
quantifiers in the neuromagnetic signals indicates the pathological state of a patient with PSE within separate areas 
of the brain. Our approach, being conveniently constructed from the set of subordinate memory functions yielding 
the rate of change of the autocorrelation function of the measured complexity dynamics, allows one to characterize 
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the neuromagnetic signals in the human brain in terms of statistical indicators. These so constructed statistical 
quantifiers in turn measure both the role and the strength of statistical memory which the underlying time series 
accommodates. Many natural phenomena are described by distributions with time scale- invariant behavior [20I ]. 
The suggested approach allows the stochastic dynamics of neuromagnetic signals in human brain to be treated in a 
probabilistic manner and to search for its statistical singularities. 

From the physical point of view the obtained results can be used as a test to identify the presence or absence of brain 
anomalies as they occur in a patient with PSE. The set of our quantifiers is uniquely associated with the emergence 
of memory effects in the chaotic behavior of the human cerebral cortex. The registration of the behavior of those 
indicators as discussed here is then of beneficial use to detect the pathological state of separate areas (sensors 10, 46, 
51, 53 and 59) in the human brain of a patient with PSE. There exist also other quantifiers of a different nature, such 
as the Lacunae's exponent, Kolmogorov-Sinai entropy, correlation dimension, etc., which are widely used in nonlinear 
dynamics and related applications, see in Ref. [21]. In the present context, we find that the employed memory 
measures are not only convenient for analysis but also ideally suited to identify anomalous brain behavior. The search 
for yet other quantifiers, and foremost, the optimization of such measures when applied to complex, discrete time 
dynamics presents a true challenge. This objective particularly holds true when attempts are made to identify and 
quantify an anomalous functioning in living systems. The present work presents such an initial step towards the 
understanding of fundamentals of physiological processes in the human brain. 

PSE is a type of reflexive epilepsy which originates mostly in visual cortex (both striate and extra-striate) but with 
high possibility towards propagating to other cortical regions [22]. Healthy brain may possibly possess an inherent 
controlling (or defensive) mechanism against this propa gati on of cortical excitations, breakdown of which makes the 
brain vulnerable to trigger epileptic seizures in patients [231 ]. However, the exact origin and dynamical nature of this 



181 ] that brain responses against chromatic 



putative defensive mechanism is not yet fully known. Earlier we showed 
flickering in healthy subjects represent strong nonlinear structures where as nonlincarity is dramatically reduced to 
minimal in patients. Here we report that patient's brain show significantly stronger statistical memory effects than 
healthy brains. A complex network composed of interacting nonlinear system with memory component is inherently 
stable and critically robust against external perturbations. Quick inhibitory effect, that is essential for the prevention 
of PSE, is made possible by the faster signal processing between distant regions. Further, such network is capable 
to facilitate flexible and spontaneous transitions between many possible configurations as opposed to being entrained 
or locked with the external perturbations 24(. In short, our findings are in line with growing body of evidence that 



physiological systems generate activity fluctuations on many temporal and spatial scales and that pathological states 
are associated with an impairment of this spatio-temporally complex structure. 
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